Newsom, J. T. (2015). Longitudinal Structural Equation Modeling: A Comprehensive Introduction. Routledge.
http://www.longitudinalsem.com/health.dat
level =~ 1* bmi1 +1* bmi2 +1* bmi3 +
1* bmi4 +1* bmi5 +1* bmi6
slope =~ 0 * bmi1 + 1 * bmi2 + 2 * bmi3 +
3 * bmi4 + 4 * bmi5 + 5 * bmi6
bmi1 ~~(vare)* bmi1
bmi2 ~~(vare)* bmi2
bmi3 ~~(vare)* bmi3
bmi4 ~~(vare)* bmi4
bmi5 ~~(vare)* bmi5
bmi6 ~~(vare)* bmi6
結構めんどい…
切片のみのモデル (model=‘no’)
線形モデル (model=‘linear’)
二次曲線モデル (model=‘quadratic’)
latent basisモデル (model = ‘latent’)
library(RAMpath)
fit.all<-ramLCM(data=data,outcome=1:6, model='all')
fit.no<-ramLCM(data=data,outcome=1:6, model='no')
fit.linear<-ramLCM(data=data,outcome=1:6, model='linear')
fit.quadratic<-ramLCM(data=data,outcome=1:6, model='quadratic')
fit.latent<-ramLCM(data=data,outcome=1:6, model='latent')
lavannコード
cat(fit.all$model$no)
## level =~ 1* X1 +1* X2 +1* X3 +1* X4 +1* X5 +1* X6
## X1 ~~(vare)* X1
## X2 ~~(vare)* X2
## X3 ~~(vare)* X3
## X4 ~~(vare)* X4
## X5 ~~(vare)* X5
## X6 ~~(vare)* X6
lavannコード
cat(fit.all$model$linear)
## level =~ 1* X1 +1* X2 +1* X3 +1* X4 +1* X5 +1* X6
## slope =~ 0 * X1 + 1 * X2 + 2 * X3 + 3 * X4 + 4 * X5 + 5 * X6
## X1 ~~(vare)* X1
## X2 ~~(vare)* X2
## X3 ~~(vare)* X3
## X4 ~~(vare)* X4
## X5 ~~(vare)* X5
## X6 ~~(vare)* X6
lavannコード
cat(fit.all$model$quadratic)
## level =~ 1* X1 +1* X2 +1* X3 +1* X4 +1* X5 +1* X6
## slope =~ 0 * X1 + 1 * X2 + 2 * X3 + 3 * X4 + 4 * X5 + 5 * X6
## quadratic =~ 0 * X1 + 1 * X2 + 4 * X3 + 9 * X4 + 16 * X5 + 25 * X6
## X1 ~~(vare)* X1
## X2 ~~(vare)* X2
## X3 ~~(vare)* X3
## X4 ~~(vare)* X4
## X5 ~~(vare)* X5
## X6 ~~(vare)* X6
lavannコード
cat(fit.all$model$latent)
## level =~ 1* X1 +1* X2 +1* X3 +1* X4 +1* X5 +1* X6
## slope =~ 0 * X1 +start( 1 )* X2 +start( 2 )* X3 +start( 3 )* X4 +start( 4 )* X5 + 5 * X6
## X1 ~~(vare)* X1
## X2 ~~(vare)* X2
## X3 ~~(vare)* X3
## X4 ~~(vare)* X4
## X5 ~~(vare)* X5
## X6 ~~(vare)* X6
fits<-round(fit.all$fit[c("chisq","df","pvalue",
"cfi","srmr","rmsea",
"aic","bic"),],
digits=2)
datatable(fits,option=list(dom='t'))
datatable(
fits, extensions = 'Buttons', options = list(
dom = 'Bfrtip',
buttons = c('copy', 'csv', 'excel', 'pdf', 'print')
)
)